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Abstract 

Based on a half-filled two-dimensional tight-binding model with nearest- 
neighbour and next nearest-neighbour hopping the effect of imperfect Fermi 
surface nesting on the Peierls instability is studied at zero temperature. Two 
dimerization patterns corresponding to a phonon vector (tt, tt) are considered. 
It is found that the Peierls instability will be suppressed with an increase of 
next nearest-neighbour hopping which characterizes the nesting deviation. 
First and second order transitions to a homogeneous state are possible. The 
competition between the two dimerized states is discussed. 
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I. INTRODUCTION 



The low dimensional electronic materials are known to be very susceptible to a Peierls 
instability towards a charge density wave (CDW) state driven by the electron-phonon in- 
teraction ffl]. The presence of a lattice distortion is usually favorable to lower the electronic 
energy and once this reduction overcomes the increase of lattice deformation energy the 
Peierls transition takes place. It has been extensively studied in a lot of quasi-one di- 
mensional (ID) materials such as organic conjugated polymers (CH) a [0 or inorganic blue 
bronzes A0.3M0O3 (A=K, Rb, Tl) as well as quasi-two dimensional (2D) materials such 
as purple bronzes AM06O17 (A=Na, K, Tl) |3|-|5| and monophosphate tungsten bronzes 
(P0 2 ) 4 (W0 3 )2 m (4 < m < 14) !§. 

It is believed that in the Peierls transition the structure of the Fermi surface (FS) plays 
an essential role. In ideal ID systems, the FS, being composed of two points seperated by 
2kp (Fermi wave vector), is always perfectly nested. The lattice distortion opens a gap at the 
Fermi level with the consequence that the energy gain from the electronic energy is always 
dominant, so that the Peierls instability with a metal- insulator transition always takes place 
(if quantum effects of phonons are not considered ||). 

The situation becomes richer in two dimensions because of a more complex FS structure. 
In general, the FS is not nested, i.e., a single mode of lattice distortion can connect only 
two points in the FS, and the gain of electronic energy from this distortion is not enough to 
overcome the increase of lattice energy. However, in some special cases the FS is still nested 
and the electronic energy may be lowered substantially by the lattice distortion even if a gap 
may be not fully opened at the FS as in one dimension. The simplest realization is the 2D 
square lattice tight-binding model with only nearest-neighbour (n.n.) hopping at half-filling, 
i.e., the 2D version of the well-known Su-Schrieffer-Heeger (SSH) model ||. In this case the 
FS consists of parallel straight lines: \k x \ + \k y \ — 7r as illustrated in Fig. 1 by the solid 
line. One of the Fermi lines may be completely moved to another by a translational vector 
Q = (tt,7t), i.e., the FS is perfectly nested with nesting vector Q. The Peierls instability for 
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this model was theoretically studied one decade ago in connection to high T c superconductors 
njnj- Unlike in one dimension, there are several possible alternation patterns for the lattice 



distortion and the corresponding bond hopping, as discussed by Tang and Hirsch ||10|| . In 
Fig. two possible dimerization patterns are shown. Both of them correspond to phonons 
with wave vector (vr,7r), which is exactly Q so that the reduced Brillouin zone boundary 
after distortion perfectly meets the original FS. The difference between them is that for case 
(a) the dimerization is in both directions, while it is only in one direction for case (b). It 
was found that even arbitrarily small lattice elastic strength will induce a lattice distortion 



into case (a) or (b), i.e, the Peierls instability is sure to occur even for this 2D model [|T0| , [TT 
On the other hand, however, it may be noticed that the present perfect nesting of the FS 
may be easily broken, for example, by introducing next nearest-neighbour (n.n.n.) hopping 
which is often not negligible. Then the following problem naturally arises: does the above 
Peierls instability still survive the imperfect nesting of the FS? 

Actually, for those quasi-2D materials which show a Peierls instability perfect nesting of 
their Fermi surfaces is never present, but an approximate, so called hidden nesting exists 
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the FS is so far away from nesting that even no hidden nesting is present. As far as real 
materials are concerned the shapes of the Fermi surfaces are obtained by band structure 
calculations and may be often rather complicated. Nevertheless, the tight-binding model 
with the n.n. and n.n.n. hopping should already be sufficient to simulate an essential 
property of real Fermi surfaces: whether they are nesting or not. Thus it is necessary to 
clarify how sensitive the Peierls instability is to the deviations of the FS from perfect nesting 
which is controlled by n.n.n. hopping. This is the topic addressed in this paper. A similar 
problem was studied previously by Lin et al. [jnj, however, their study was only limited to 
the dimerization pattern (a). (Actually this pattern is not favorable in a large region of t' 
as will be seen later.) In addition, the main result Fig. 7 in their work was not convincingly 
presented. We will include the two possible patterns (a) and (b) and address the unexpected 
competition between them. 
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FIG. 1. The Fermi surfaces for the 2D tight-binding model without dimerization at half-filling. 
The solid line has if = 0, the dashed line has t' = 0.1, and the dotted line has t' = 0.2. 



FIG. 2. The lattice distortion patterns (a) and (b). In the figure a thick solid line corresponds 
to a strong bond with hopping integral t(l + 5), a dashed line corresponds to a weak bond with 
hopping integral t(l — 5), and a thin solid line corresponds to a normal bond with hopping integral 
t. Both patterns correspond to phonons with wave vector (tt,tt). The dimerization is along two 
axes for case (a), while only along x axis for case (b). 



We begin with the following Hamiltonian based on a square lattice with a half-filled 
band: 



(a) 



(b) 



II. FS NESTING AND PEIERLS INSTABILITY 
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H = ~t 121 1 + a ( U tj - U i+l,M C h,<rCi+l,j,a + h-C.) -tJ2i 1 + a ( U lj ~ U lj + M C lj,aCi,j + l,a + h.C.) 

K 

12 i C l,j,aCi+l,j+l,a + cJj, ff Ci+lJ-l, ff + h - C + y li''"., - Wf+lj) 2 + «j ~ <j+l) 2 ] > (!) 

where cj J)0 .(cj JiCr ) denotes the creation (annihilation) operator for an electron at site 
with spin a (i denotes x coordinate and j denotes y coordinate), u^j V represents the displace- 
ment component of site in x/y direction, t, t' are n.n and n.n.n. hopping parameters 
and a is the electron-lattice coupling constant. The last term above describes the lattice 
elastic potential energy with K the elastic constant. The lattice kinetic energy is omitted 
here since we do not study the dynamic behavior of phonons. For the lattice distortion 
patterns investigated here, the distance between n.n.n. sites remains unchanged, therefore 
no dimerization of t' is considered in Hamiltonian (O. 

The lattice distortion, as shown in Fig. 2 having the wave vector Q, may be explicitly 
introduced as 

<, - »•;•.,., = (•-! )'•'"• "'!, - "tj . , = i ir-'" 

for case (a) and 

<, - »;.„ = ( ir- ? "- <, - "1, . , = o 

for case (b), where u is the amplitude of dimerization determined by minimization of the 
ground state energy. For convenience, two dimensionless parameters are defined as follows: 
the dimerization amplitude 5 = au and the electron-lattice coupling constant rj = a 2 t/K. 
The n.n. hopping integral t is taken as the energy unit. 

Before proceeding, we would like to mention that two further possible lattice distortion 
patterns as discussed in previous works on the 2D Peierls instability, which correspond to 
wave vector (vr, 0) and/or (0,7r) [fTCO , are excluded in our study. We omit them because 
at t' — their energy gains are much smaller than those for the patterns considered here 
|TTJ . (They are possibly favorable only when a large Hubbard U term is switched on [|HJ.) 
This is physically understandable: they have a different wave vector from the nesting one. 
Very recently, a more complex lattice distortion pattern with incommensurate structure was 
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studied by Ono and Hamano in the case of t' — |L4| . However it will not be included here 
since it is not unique [0. The patterns considered in Fig. 2 may be regarded as typical 
structures for the study of the Peierls instability in two dimensions. In the following the 
Peierls instabilities for t' = and t' ^ are discussed respectively. 

A. t' = 

For completeness, we first reproduce the results for the perfect nesting case t' = 0, which 
are straightforward but helpful. In momentum space the electronic spectra for case (a) and 
(b) are respectively written as 



a = ±2\/ (cos k x + cos k y ) 2 + <5 2 (sin k x + sin k y ) 2 , 



e k b = ±2y (cos k x + cos k y ) 2 + S 2 sin k x . (2) 

In the ground state only the lower band is fully occupied for each case (i.e., the chemical 
potential is zero). Then the ground state energy is given by 

E = 2^2e Ka + N5 2 / V 

k 

for case (a) and 

E = 2Y,s Kb + N5 2 /2 V 

k 

for case (b), where the summation is over the wave vector k = (k x , k y ) in the Brillouin zone 
—7i < k x ± k y < it and N is the total number of lattice sites. The factor 2 in front of the 
summation is due to spin degeneracy. Obviously the increase of elastic energy for case (a) 
is two times of that for case (b). 
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FIG. 3. The optimal value 5* as a function of rj for t' = 0. The dashed line is for case (a) and 
the solid line is for case (b). 
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FIG. 4. The energy gain of the dimerized state: E* at 5 = 5* minus Eq at 5 = as a function 
of rj, corresponding to Fig. 0. 



The results are shown in Figs. |3] and |j. In Fig. |3] the optimal dimerization parameter 
5* with lowest energy vs. the electron-lattice coupling 77 is plotted, and Fig. f| gives the 
corresponding ground state energies E* at these S* values (with respect to the energy of 
the undimerized lattice Eq which is equal to —16/n 2 times N). From these figures it can 
be seen that the Peierls instability takes place as long as the electron-lattice coupling is 
non-zero for either case (a) or (b). Although the magnitude of 5* is small when 77 — ^ 0, it 
is proven to be finite, with an exponential dependence of exp(— c/rj) on 77 (c is a constant). 
This confirms the related statement made in the Introduction. Moreover, it is seen from 
Fig. § that the energy E* for case (b) (solid line) is always lower than the corresponding one 
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for ) (dashed line) in the full region of finite 77, which means that the case (b) is the 

more favorable dimerization pattern for if = 0. Unfortunately in their works, Mazumdar as 
well as Tang and Hirsch incorrectly prefered case (a), see Ref. A similar conclusion as 
ours was reached by Ono and Hamano [|l4|j . 

B.^0 

For if 7^ the perfect nesting of the FS is broken. In Fig. [I] the Fermi surfaces for the 
undimerized lattice are plotted for several if values. One may see that the bigger t' is the 
farther the FS deviates from perfect nesting. 

Still the Hamiltonian can be easily diagonalized, but the final results are quite nontrivial. 
Now the electronic spectra become 

£ k,a = ~ 4t' cos k x cos k y ± 2 J (cos k x + cos k y ) 2 + <5 2 (sin k x + sin k y ) 2 , 

£^,6 — — 4t' cos k x cos ky ± 2 J (cos k x + cos k y ) 2 + 5 2 sin 2 k x , (3) 

and the ground state energy is 

E = 2 ]T e± a + N5 2 h 

for case (a) and 

E = 2 ]T e± b + N5 2 /2 V 

for case (b), where fi a and fib are chemical potentials for ) and (b), respectively. By 

solving the equation BE/ '85 = we may single out the optimal value 5* for each t' with 
fixed 77. In the following a typical value rj = 0.5 is adopted |T3|JT^| . 

The results for 5* are shown in Fig. ^. It is clearly seen that 5* goes to zero with 
increase of t' in both cases, which means that the Peierls instability is suppressed at some 
t', as expected. And more interestingly, the details for both cases are different. For case 
(a) 5* first decreases weakly with increase of t', but at some critical t' cl ~ 0.1733, it drops 
suddenly to zero, showing a first-order transition. On the other hand, for case (b) 5* first 
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retains its t' = value, and after t' is beyond about 0.12 it begins to decrease gradually and 
approaches zero smoothly at t' c2 ~ 0.1704 — the transition is of second-order. The value 
t' c 2 is close to, but different from t' c l . 

0.25 

0.2 
to 0.15 

0.1 
0.05 

0.05 0.1 0.15 0.2 
t' It 

FIG. 5. The optimal value 5* as a function of t' for r] = 0.5. The dashed line is for case (a) and 
the solid line is for case (b). 

To see the above transitions more clearly, we plot in Figs. || and [7| the dependence of the 
ground state energy on the dimerization parameter 5 for several different t' values for case 
(a) and (b), respectively. For ) it is seen that when t' approaches the critical value 

t' cl , the energy E first increases with 5 and then decreases towards a minimum, as shown 
in the inset of Fig. [| Thus this minimum may be only a local one. It is expected that 
at t' cl the energy E at such a minimum is the same as that at 5 = 0, see the inset. And 
then once t' is beyond t' c x the 5 value with absolutely lowest energy should be taken zero, 
i.e., 5* = 0. So a first-order transition arises. On the other hand, for case (b) the energy 
E always decreases first with increase of 5 until a minimum is reached. So this minimum is 
actually a global one. It shifts continuously towards zero with increase of t', which explains 
the second-order transition. 
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FIG. 6. The ground state energy per site as a function of 5 for several t' values for case (a). 
The solid curves from up to down correspond to t 1 = 0.164, 0.168, 0.172, 0.1733, respectively, and 
the dotted line connects the global minima of them. The curves for t 1 = 0.172 and t' = 0.1733 are 
enlarged in the inset, where the initial increase in each curve is shown. 
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FIG. 7. Same as Fig. || for case (b). The curves from up to down correspond to 
t' = 0.14, 0.15, 0.16, 0.17, respectively. In the inset the curve for t' = 0.17 is shown in the 
small 5 region. 

It deserves to point out the physical reason for the constant value of 5* in the small t' 
region (about t' < 0.12) for case (b). Actually, in this region the two energy bands e^ b do 
not overlap when one takes relatively large 5 value like 5*, so that the lower band e^ b is fully 
occupied and the upper band e^ b is empty in the ground state. Then the electronic energy 
is simply given by Sk^k^ which is independent of t' due to J2k cos k x cos k y = 0. Therefore 
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the ground state energy for t' 7^ is the same as that for t' = 0, so is the solution S*. 

While so far the results of the Peierls instabilities for the two patterns have been given 
separately, it is now natural to think of the competition between them. The competition 
becomes evident when the respective lowest energies E* at different t' values are compared, 
as shown in Fig. [5[ It is interesting to see that the energy E* for case (b) is first less than 
that for ), but later becomes larger with increase of t' . Consequently, a transition 

between these two dimerization patterns is predicted. To summarize, the evolution of the 
stable state for the system with increase of t' may be described as follows: at t' = the 
dimerized state stabilizes due to the perfect FS nesting, and moreover the dimerization 
pattern takes form (b). This state remains stable until at some t' value (which is about 0.13 
for 7] = 0.5) it is replaced by the other dimerized state (a). Finally, the dimerized state 
breaks down, i.e., the Peierls instability is suppressed at some critical value t' c which is equal 
to t' c j as can be seen in Fig. |8|. Thus the suppression is a first-order transition. 

-1.628 
-1. 63 
^- -1. 632 
-1.634 
-1.636 

0.05 0.1 0.15 0.2 
t' It 

FIG. 8. The energy at 5 = 5* as a function of t' for rj = 0.5, corresponding to Fig. ||. Both 
curves stop at their respective critical points. 

Above, only the results for t' > are presented. What is expected for t' < 0? Actually, 
taking t' — ► —t' reflects all the results as shown by Figs. f| and § about t' = 0. This 
conclusion may be proven as follows. From Eq. (0) one may get the relations between the 
electronic spectra and their corresponding ones (e^) after reflection t' — > —t': = —e^. 
(Because of the irrelevance of the proofs to the dimerization pattern the subscript which 
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differentiates case (a) or (b) is omitted.) Then one obtains that the chemical potential p, for 
the spectra is just the negative of that for the e^, i.e., JX = —(jl at half-filling, due to the 
following equations: 

E 1= E 1= E 1= E i = #/2. 



-4<-f 



In addition, one needs to notice that the summation for all energy levels Dk^k 
X)k ~~ 8t' cos fca, cos fcy = 0. Finally the following relations hold: 



E< 



E^ 



E 



E ( 



Note that the first and the last terms in the above formulas are nothing but the electronic 
parts of the ground state energies for t' and —t', respectively. They are proven to be the 
same, so are all the subsequent results derived from them. 

Up to now all results have been shown for fixed rj. In the last part of this section we 
will look at the role of rj. Since the parameter 77 only appears in the elastic energy, i.e., it 
is irrelevant to the electronic part of the ground state energy, the qualitative properties of 
curves E vs. 5 as shown in Figs. [| and |7| are expected to be retained, so are the above 
qualitative results for the Peierls instability. Quantitatively one may think, the smaller the 
electron-lattice coupling rj is (or equivalently the larger the elastic strength K is), the smaller 
the critical value t' c for the suppression of the Peierls instability should be. As a check, the 
case for rj = 0.3 is calculated. The results are shown in Figs. |9| and |10|, which confirm the 
above predictions. 




12 



FIG. 9. The same as Fig. | but with rj = 0.3. 
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FIG. 10. The same as Fig. [8| but with rj = 0.3. 



III. DISCUSSION 



In the previous section the Peierls instability for two dimerization patterns has been 
carefully studied in its dependence on t'. The results are instructive. It is now known that 
low-dimensional metals show two types of electronic instabilities: either a Peierls instability 
or a superconducting one. The Peierls instability may prevent some metals from entering a 
superconducting state. For example, in the transition metal bronzes, the Peierls instability is 
the dominant mechanism and superconductivity was only found in LiMoeOi7 Our results 
suggest that the Peierls instability may be suppressed by some way of increasing imperfect 
nesting, for example, by applying pressure to enhance the n.n.n. hopping. Moreover, the 
possible transition between different dimerized states is reminiscent of the experimental 
results in monophosphate tungsten bronzes like P4W12O44, where double Peierls transitions 
occur with change of temperature |||7j . The competition between two dimerization patterns 
at finite temperatures is believed to be interesting and will be left for future investigation. 

A further problem is the consideration of effects of electron correlations, e.g., the Hubbard 
U term which is not included in this work. For t' — the problem has been studied by 



Tang and Hirsch by numerical calculations [10,11 . It was found that the on-site Coulomb 



13 



interaction weakens the dimerization in two dimensions as soon as U is present [fn|| . The 
possible explanation is that the U term favors the appearance of antiferromagnetic (AF) spin 
order in 2D half-filled Hubbard model, while the dimerization stabilizes local spin singlets 
which are unfavorable for strong U. Above, we showed that t' also suppresses dimerization. 
Does this indicate that the simultaneous presence of both U and t' will speed the suppression 
of the dimerization? We leave this issue for future work and only give some clues here from 
previous analysis in the large U limit (for the moment the concrete dimerization pattern is 
not considered). In this limit we are facing a problem of localized electrons interacting via 
an effective n.n. exchange J (~ t 2 /U) and n.n.n. exchange J' (~ t' 2 /U). The Peierls system 
then transforms into the corresponding spin-Peierls (SP) system. The 2D (or quasi-lD) 



spin-Peierls instability without n.n.n. exchange J' was studied by some authors fI5|Jlq| . It 
was found that the SP transition does not spontaneously occur unless the so called spin- 
lattice coupling (analogous to 77 here) exceeds a threshold. On the other hand, it was known 
that the n.n.n. exchange J' will frustrate the AF order for the 2D spin system. So one may 
assume that the additional inclusion of J' will be favorable to the formation of a SP state. 
Thus it is not unreasonable to believe that the effects of t' and U on the Peierls instability 
may cancel in part when acting simultaneously, although each of them separately tends to 
suppress it. 



IV. CONCLUSION 

In conclusion, a 2D tight-binding model with n.n. hopping t and n.n.n. hopping t' is 
used to study the effect of imperfect FS nesting on the Peierls instability of the ground state. 
Two possible dimerization patterns corresponding to a phonon vector (tt, tt) are considered 
as case (a) and (b). It is found that the Peierls instability will be suppressed with an increase 
of t' which characterizes the deviation from perfect nesting. The details for the two cases 
are different: for case (a) the suppression is a first-order transition while for case (b) it is of 
second-order. Also a transition between the two dimerized states is investigated. 
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